Bhabha scattering at NNLO with next-to-soft stabilisation

A critical subject in fully differential QED calculations originates from numerical instabilities due to small fermion masses that act as regulators of collinear singularities. At next-to-next-to-leading order (NNLO) a major challenge is therefore find stable implementation numerically delicate real-virtual matrix elements. In the case Bhabha scattering this has so far prevented development fixed-order Monte Carlo at NNLO accuracy. paper we present new method for stabilising element. It based on expansion soft photon energies including non-universal subleading term calculated with regions. We have applied obtain and efficient within McMule framework. first time results photonic corrections scattering.